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ZEROS OF RANDOM FUNCTIONS GENERATED WITH 
DE BRANGES KERNELS 

JORGE ANTEZANA, JORDI MARZO, AND JAN-FREDRIK OLSEN 


Abstract. We study the point process given by the set of real zeros of 
random sums of orthonormal bases of reproducing kernels of de Branges 
spaces. Examples of these kernels are the cardinal sine, Airy and Bessel 
kernels. We find an explicit formula for the first intensity function in 
terms of the phase of the Hermite-Biehler function. We prove that the 
first intensity of the point process completely characterizes the underly¬ 
ing de Branges space. This result is a real version of the so called Calabi 
rigidity for GAFs proved by M. Sodin. 


1. Introduction 


1.1. Background. Let E(z ) be a function of the Hermite-Biehler class, i.e., 
E(z) is entire and satisfies the inequality |E(, 2 )| > \E*(z)\ for Im z > 0, 
where E*(z) = E(z). Given such a function, the de Branges space H{E) is 
defined by 

H ( E ) = {/ entir e : G # 2 (C+) j , 

with norm given by 



Here, we use H 2 ( C+) to denote the Hardy space of the upper half-plane. 
The standard reference for de Branges spaces is the book [9] by de Branges. 

Note that in this paper, we restrict ourselves to considering Hermite- 
Biehler functions without real zeroes. This excludes the existence of points 
i£R such that f(x ) = 0 for all / E H(E). 

The spaces H(E) can be thought of as weighted versions of the Paley- 
Wiener spaces PW : 2 . The simplest description of PW 2 is perhaps as con¬ 
sisting of the Fourier transforms of functions in L 2 (— a, a) with the induced 
norm of L 2 (M). By the Paley-Wiener theorem, this space is identical to H{E) 
with E{z) = e -m2 . And, as we point out below, its reproducing kernel is 
given by translates of the cardinal sine function, i.e., 


Ky(x) 


sina(x — y) 
{x - y ) 


Other important examples are obtained by using the Airy or Bessel func¬ 
tion. The corresponding reproducing kernels, together with the cardinal sine, 
are ubiquitous in random matrix theory. A reason is that the determinantal 
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random point processes given by these kernels (see below for a definition) 
describe the possible asymptotic distribution of eigenvalues from large Her- 
mitian matrices, see, for instance, [7J. 

In this paper, instead of considering determinantal processes, we investi¬ 
gate random point processes whose relation to the cardinal sine, Bessel and 
Airy functions is slightly different. Namely, translates of these functions es¬ 
sentially appear as orthonormal bases of reproducing kernels of certain de 
Branges spaces. We study the point processes obtained by considering the 
random linear combinations obtained from each of these bases. These are 
examples of Gaussian Analytic Functions (GAFs). 

Remark 1.1. The de Branges spaces can be seen as a subclass of model 
subspaces of the Hardy space H 2 ( C+). The model subspaces Ad(@) of H 2 
are exactly those of the form _Ad(0) = H 2 n (0H 2 ) 1 , where 0 is an in¬ 
ner function. That is, 0 is bounded on the upper half-plane and satisfies 
|©| = 1 on M. To see the connection, one may easily verify that if E is 
of Hermite-Biehler class, then © = E* /E is a meromorphic inner function, 
and, conversely, it is possible to show that any meromorphic inner function 
can be factorized in this way. In other words, the operator / i —> f/E, where 
0 = E*/E, is an isometry from H{E) onto _M(0). See [12] for more details. 


To introduce the GAFs that are our object of study, we first recall some 
additional facts and notations about the spaces H(E). For a point w E C, 
the reproducing kernel function K w E H(E) is the function satisfying 

f+°° _ dx 

f(w) = {f , K w ) = J f(x)K w (x) | 2 

for all / E H(E). This function exists since point evaluations are bounded 
functionals. A straight-forward computation yields the formula 


K w (z) = K(z,w) 


i E{z)E(w) - E*(z)E*(w) 
2n z — w 


It is useful to introduce the polar decomposition E(x) = \E(x) [e -1 ^^ for 
x E M. Here, the so called phase function (ft is an increasing C°°(M) function. 
With this notation, for x E M, we have 

II AT X || 2 = K(x,x) = -fi(x)\E(x)\ 2 , 

7 r 

and, for x, y E M, we get the following expression for the normalized repro¬ 
ducing kernels 

K(x,y) _ \E(x)\ sin(</>(x) — (f>(y)) 

V[X) K(y,y )V2 x - y 

Observe that k y (x ) = (k y , k x ) = 0 whenever x, y E M are such that </>(x) — 
4>(y) = kir for fc E Z. This means that if {w n } is the sequence of points such 
that <j>(uJ n ) = a + 7rn, for some a E [0, n), then the family of functions {k Un \ 
forms an orthonormal system in H(E). However, as de Branges proved, see 
m p. 55], more is true. 
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Lemma 1.2. Suppose that H(E ) is a de Branges space. Then, for all 
a E [0,7r), except at most one, the system {k Un } is an orthonormal ba¬ 
sis for H(E). Moreover, these are the only orthonormal bases of repro¬ 
ducing kernels, and the exceptional a is characterized by the condition that 
e iaz E(z ) - e~ iaz E*(z) E H(E). 

We now give the definition of a de Branges GAF. 


Definition 1.3. Let H{E ) be a de Branges space with phase function <j>(x) 
and an orthonormal basis of reproducing kernels {k Un }. For real i.i.d. standard 
normal random variables a n , we call the random function 


F ( x ) =^2 a nkoj n (x) = 

n n 

-E 


i E(x)E(u n ) - E*(x)E*(u n ) 
"2 yjTuf'[x)\E{x)\{x - u n ) 
\E(x)\ sin - (j>{u n )) 

&n / 

y/irft (uj n ) x - u n 


a de Branges GAF. 


(1) 


We remark that Kolmogorov’s inequality shows that, for a fixed x E M, the 
series (jlj) is almost surely pointwise convergent. In fact, since Y2 n \k Un (z)\ 2 
converges uniformly for z on compact subsets of C, it can be proved that 
the series in ([!]) almost surely converges uniformly on compact subsets of C. 
Thus, the de Branges GAF is an entire function with probability one. 

More generally, for a sequence of analytic function {'tp n (z)} defined in a 
certain region A C C, if ^ \ip n (z)\ 2 converges uniformly on compact sets of 
A and a n are real i.i.d. standard random variables, the series 

^(z) = YanM*) ( 2 ) 

n 

defines an analytic function with probability one. If a n are complex i.i.d. standard 
normal random variables, the function is called a Gaussian Analytic 
Function (GAF), see [13| Lemma 2.2.3]. 

The main feature connecting the zeros of the de Branges GAF F(x) with 
the de Branges space H(E) is that, by definition, the covariance kernel of the 
gaussian stochastic process (F(x)) x coincides with the reproducing kernel 
of the de Branges space H(E ), i.e., 

E [F(x)F(y)\ = K(x,y). 

Here E denotes the expectation. Recall that the covariance kernel, as (F(x)) x£ ^ 
has mean zero, completely defines the behavior of the process. 


Remark 1.4. By Remark 0 above, the random function F/E should be 
called a model space GAF. Since E{x) is a non-random function without 
zeros on the real line, the zeros of F/E has exactly the same distribution as 
that of F. We will therefore pass freely between the two. Note that if we 
factor out non-random factors, then we can write 


F(x) 

E(x) 


2 1( ^(*) sin — a) 


7r 


E' 


(-i)’ 


Y^KJ^-Wn)’ 


where a E M is the parameter so that {w n } is the sequence of points such 
that 4>(u> n ) = a + 7rn. 
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To formulate our results, we need to define the intensity functions, see 
na Chap. 1], which describe the distribution of a point process and, in 
some sense, are similar to densities. To this end, suppose that a simple point 
process X is defined in M or C, i.e., the points of X are almost surely of 
multiplicity one. In our case, this process will be the real zeroes of a de 
Branges GAF ([!]). The joint intensities p k (xi ,..., x k ) are defined by the 
relations 

E[X(D 1 ) ■ ■ ■ X(D k )\ = / p k (xi,...,x k )dxi...dx k , 

J Dix---xD k 

where X(A) stands for the number points of X in a Borel set A C M, the 
Borel subsets D i,..., D k C M are mutually disjoint, and p k {x i,.. . , x k ) = 
0 if Xi = Xj for i ^ j. Other reference measures can be considered, but 
throughout this work we will consider only intensity functions with respect 
to the Lebesgue measure. 

In particular, if Z^(F) denotes the set of real zeroes of a de Branges GAF 

F (x) = ^2a n k Un (x), 

n 

its first intensity satisfies 

E[#(Zr(F) n D)\ = f Pl (x)dx. 

Jd 

A particularly well studied class of point processes is defined by assuming 
that the intensities are given by determinants of certain kernel functions 


Pk(x i, ...,x k ) = det(K(xi,Xj))ij = i t _ jk . 


These are exactly the deteminantal point processes mentioned above, see 
[2 EE]. The de Branges reproducing kernels provide some well known ex¬ 
amples of determinantal point processes. These processes have been studied 
specially in connection to random matrix theory. More specifically, the car¬ 
dinal sine kernel describes the distribution of neighbouring eigenvalues in the 
bulk for the Gaussian unitary ensemble, GUE [20], the Airy kernel (13) de¬ 
scribes the soft edge of the spectrum for the GUE [30], while the Bessel kernel 


(18) represents the hard edge of the spectrum of the Jacobi and Laguerre 
ensembles, [31]. 


In the determinantal case, the first intensity function is just given by the 
diagonal of the kernel, namely pi(x) = K(x,x). For GAFs defined as in 
<© with the a n being complex i.i.d. standard normal variables, Sodin |26| . 
computed that 

pi (z) = -A log K (z,z) (3) 


7T 


whence 

E[#(Z C (V) n D)\ = - [ Alog K(z,z)dz 

7T J d 

where D C A and Zc( v k) is the zero set of T^). This generalises previous 
works of Kac |15| . Rice HQ [25], and Edelman-Kostlan [10] , 

In our case, as we consider real random coefficients, and not complex, the 
first intensity cannot be computed using the expression ([3|. Instead, the 
situation is more similar to the case of real zeroes of random polynomials 
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with real coefficients, studied in |29j . Indeed, in the case of real random 
coefficients, the de Branges GAFs have zeroes which are symmetric with 
respect to the real line (recall that the kernel functions are symmetric with 
respect to the real line), and one expects to find a certain proportion of 
them on the real line. For this reason, with respect to the complex plane, 
the intensity functions should have a singular part that is supported on the 
real line. Indeed, this is a consequence of the following formula which was 
proved in HD. It says that in the distributional sense, the first intensity 
function satisfies 

Pi(z) = ^ A1 °g (k(z,z) + y / K(z,z) 2 - \K(z,z)\ 2 ^J . (4) 

1.2. Results. In the following theorem, we give a formula for the first in¬ 
tensity function of de Branges GAFs in terms of the Schwarzian derivative 
of the phase function <f>, i.e., 



Theorem 1.5. Let H(E ) be a de Branges space with phase function cf. For 
any choice of orthonormal basis of reproducing kernels, then the first intensity 
function of the corresponding de Branges GAF satisfies 


Pi 0*0 


1 / (//(x) 2 

vr V 3 ' 


S <j>{x) 
6 


To compute p\(x) one can use any of the following equivalent formulas: 


Pl (x)=E[\F'(x)\ : F{x) = 0] 

1 / S 2 

= ~\ wtw- logAI(s,t) 

7T V OtOS s=t=x 

_ sup {| h! (x)| : \\h\\ H (E) = h( x ) = °} 
7 Ty/K(x, X ) 


(Rice formula) 
(Edelman-Kostlan) 
(Bergman metric), 


where E[|T'(x)| : F(x) = 0] denotes a conditional expectation. 

We found the appearance of the Schwarzian derivative, which is an in¬ 
finitesimal version of the coss-ratio (see w quite intriguing. It is invariant 
under any Moebius transformation and it measures, in some sense, how far 
a map is from being a Moebius map. Geometrically, it is connected with 
curvature. It appears in the study of univalent functions, generalizations of 
Schwarz-Christoffel mappings, Sturm-Liouville equations and real dynamics, 
see the survey paper [22] . 


Example 1.6. The Paley-Weiner space PWf is a de Branges space with 
E(z) = e~ mz . It follows that <f(x) = ax, and that the Paley-Wiener GAF is 
given by 


F (x) = 

n£Z 


tin 


sina(x — n) 
a(x — n ) 


In t his case, the process given by the real zeroes is stationary. By Theorem 
1.5 it follows that pi(x) = a/iry/3. Although well-known, we point out 


that this is surprising since it implies that the expected number of zeroes 
to be found in an interval is strictly smaller than what one would expect 
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by comparing to the density of the zeroes of the reproducing kernel. Some 
properties of this process where studied in [4J. 

Below, in Section [3j we will consider examples of de Branges spaces given 
by the Airy and Bessel functions. Here, we state the following lemma which 
will be used to transfer information of the derivative of the phase function 
4>'(x) to the first intensity function pi(x). 

Lemma 1.7. Suppose that 4>(x) is a function such that 

(j 

where g(x) = o(x~ a ), g'(x) = o(x~ a ~ 1 ) and g"{x) = o(x~ a ~ 2 ) as x -A oo, 
respectively. Then, as x —>• oo, 

Sm = ^pi + o^). 

Proof. The proof follows from a direct computation using the definition of 
Sg. m 

An interesting feature of GAFs is that the distribution of the zero set 
depends just on the first intensity function, as was shown by Sodin in [26] 
(see also |13| Section 2.5.]). The idea is that from ([3|, it can be shown that 
the difference between the logarithms of the covariance kernel of two GAFs 
with the same first intensity is an harmonic function. This implies that the 
covariance kernels are equal on the complex diagonal, up to a multiplication 
by a non-vanishing harmonic function. From this, one is able to argue that 
the kernel functions are essentially equal even off the diagonals, which leads 
to the desired result. 

As Sodin mentioned in his paper, his result can be seen as a special case of 
Calabi’s rigidity theorem, [8]|I6|. A version in the setting of model subspaces 
was proved by Nikolskii, [21] , In our case, as we are working in a real context, 
the formula ([3]) is no longer valid. However, as we show in the following 
theorem, the first intensity still determines the distribution. 

Theorem 1.8. Let F(x), G{x) be two de Branges GAFs and let pf (x), pf (x) 
be the respective first intensity functions. If, for all x £ M, 

Pi( x ) = Pi ( x ), 

then there exists a non-random, analytic function S(z), which does not van¬ 
ish anywhere, such that S(x)F(x) and G(x) have the same distribution. In 
particular, the random processes given by the zeros of F(x) and G(x) have 
the same distribution. 

This result follows as a consequence of the following theorem, which shows 
that the first intensity function not only determines the de Branges GAF, 
but also the underlying de Branges space up to a very special isometric 
isomorphism. 

Theorem 1.9. Let E\(z), E 2 (z) be Hermite-Biehler functions with phase 
functions respectively given by (f>i(x), (feix) and let K\(z, w), K 2 (z, w) repro¬ 
ducing kernels of the de Branges spaces H{E\), H(E 2 ). If 

pf x (x) = pf 2 (x), for all x £ R, 
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then there exists an entire function S(z ) without zeros, such that 
K 2 (z,w) = S(z)K 1 (z,w)S(w), 
and the map f i—>• Sf is an isometry from H(E\) to H(E 2 ), 

Observe that it follows from this result that a de Branges GAF is not 
stationary unless the kernel is essentially the reproducing kernel of a Paley- 
Wiener space, (see (T3J Corollary 2.5.4.]). 


1.3. Structure of the paper. In Section 2, we prove Theorem 1.5 In 


Section 3, we apply Theorem 1.5 to several examples, including de Branges 
spaces with the cardinal sine, Bessel and Airy functions yielding the re¬ 
producing kernel as well as to de Branges spaces connected to orthogonal 
polynomials. In Section 4, we prove Theorem |1.8| and Theorem |1.9| 


1.4. Notation. We use the notation f(x) < g(x ) to indicate that f(x)/g(x) 
is bounded above by some positive constant. We write f(x) ~ g(x) if both 
f{x) < g{x) and g(x) < f(x) hold. 


2. The first Intensity function 


In this section, we explain two methods for computing the first intensity 

The first method 


function of a de Branges GAF. This proves Theorem 1.5 


relies more explicitly on probability theory, and is based on the Rice formula 
(see formula ([5]) below). The second method is function theoretic and is 
based on solving a deterministic maximal problem. This maximal problem 
is arrived at via the classical Edelman-Kostlan formula and its connection 
to the Bergman metric. 


2.1. Computation of pi{x) using Rice’s formula. As we consider Hermite- 
Biehler functions E(z) without zeros on M, the process given by the zeros of 
the de Branges GAF F(x), defined by Q, has the same distribution as the 
process given by the zeros of 


F 0 (x) 


F ( x ) = i</,(*) y- sin(0(a) - <j>(u n )) 

E ( X ) n U - UJ n ) 


^ ' ®ri^0 {x , Wn)i 
n 


where ko(x,uj n ) is the normalized reproducing kernel in _A4(0) (see Remark 
0 >- Observe that this corresponds to considering the GAF defined by an 
orthonormal basis of reproducing kernels of the model subspace A4(0). 

According to the Rice formula (see [2] Chap. 11]), the first intensity of 
the point process associated to the real zeroes process Fq(x ) is given by 


pi(x) = E [|F4(x)| : F e {x) = 0] = [ \t\p F {x)iF m x) (0,t)dt, (5) 

Jm. u 

where P F q( x ),F' 0 (x) is th e joint probability density of the two dimensional 
normal vector ( F e (x ) X (x)). Since both F e (x) and F, 4 (x) are normally 
distributed, it follows that 

PF e (x),F^x)(s,t) = 27r ^pi e xp(^(8,t)S- 1 (8,t) T ), 
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where the covariance matrix £ is given by 

( E[F e (x) 2 ] E[F e (x)i%(x)]\ fa b 


E[F e (x)F4(x)] E[F4(.x) 2 ] 


c d 


Let Kq[x,o) be the reproducing kernel in _M(0). We compute the entries: 
a = E [. F e (x ) 2 ] = y,ke(x fUn ) 2 


= Kq(x, x) = 


O c 

d = E [. Fq(x ) 2 ] = ^ (^—fc©(x,w n )^) 

n 

d 2 

= LrNw A "e(s,i) 

OSOt s=t=x 


2 (j)'(x) 3 + <j)"'(x) 
6n 


6 = c = E [F 0 (x)F4(x)] = ^ ke(x,u n )—ke(x,ujn) 

n 

n 

“2 9x e{ ’ 2vr ' 


On the other hand, since 


it follows that 


= - 6 


|£| y - c a 


Pi(*) = 


27TISI 1 / 2 J 

|£|V 2 


| e -l(0, t )E-W dt 


00 at 2 

te 2 i s i dt 


As | £ | = ad — be, this yields 

pi(x) = -J-(d- — 
it V a V a / 


1/7T / 2(0'(x)) 3 + (j)'"(x) TT (</>"(x)) 2 


7T y </>'(x) 


1 4>'(x) 2 Scj)(x ) 

7T V 3 6 


^'(x) 47T 2 


where, as before, denotes the Schwarzian derivative of (j). 
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2.2. Computation of p\(x) using the Bergman metric. To compute 
the first intensity function, one can alternatively use the Edelman-Kostlan 
formula, |10| . It says that for a real point process with covariance kernel 
K(s,t), it holds that 


Pi(x) = 


1 


d 2 


7 r V dtds 


log K(t,s) 


t=S=X 


( 6 ) 


Following Bergman [6] p. 35, formula (27)], one can relate this expression 
with an extremal problem 

d 2 . , inf {ll^ll^(E) : K x ) = l} 2 


dtds 


log K(t,s) 


s=t=x K ^ X,X ^ [ni {\\ h \\H(E) '■ h ( X ) = °> h '( x ) = !}’ 


where it is implicit that both infinrums are taken only over h £ H(E). This 
can be reformulated as 


d 2 

dtds 


log K(t,s) 


S=t=X 


sup{|h'(x)l 2 : ; /t ]//(/■;) = 0, h(x) = 1} 
K(x, x ) 


(7) 


This last expression is related with the Bergman metric for the space H(E) 
(see e.g., m Ch. 6] and |5]). Indeed, the first intensity p\{x) is the density, 
with respect to the real line, of the Bergman metric given by the reproducing 
kernel of H(E). Using ([7]), we can compute p\(x) by solving this (determinis¬ 
tic) extremal problem in H{E). To this end, we prove the following technical 
lemma. Just as a curiosity, this lemma provides a generalisation of the Basel 
problem. 


Lemma 2.1. Let H(E) be a de Branges space with phase function (j). If 
{kui n } is an orthonormal basis of reproducing kernels for H(E), then 

^ 4>'{uk) = ^Vfc) 2 + Sct>{u k ) ^ 


n^k 


(u k - w n ) 2 ()'(w n ) 


6 


Proof. Fix n, and let x G M \ {ui n }. Computing the norm of the normalized 
reproducing kernel k x . we get 

tf(x) 


ET7 


1 


whence 


M 


(j>'(u k ) 


4>'(u n )(x - ui n ) 2 sin 2 (a - <f(x)) ’ 


<j>'{u k )(f)'{x) 


n^k 


(j)'(u n )(x - u n ) 2 sm 2 (<f(u} k ) - (j)(x)) (x - Uj k ) 2 


ft {ukW (x) 


tf{ukW{x) 


+ 


,sin 2 (<^( u k )-(j>(x)) (4>{oj k ) - 4>(x)) 2 

f (j>'(uj k )(f)'(x) _ 1 

\((/>(u k ) - 4>(x)) 2 ( x-u k ) 2 

The first summand converges to (j)'(uj k ) 2 /3 as x — > ix k since 

«m( 1 -V 1 

x->o \sin (x) x z 
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To deal with the second summand, we use the well-known formula 


lim 

x,y-H 


( 4f{x)4f{x) 
VOO) - <Xy)) 2 


(x - yf 


\sm- 


It follows from Lemma |1.2| that for almost every x E M, there exists a 
sequence {w n } C M such that {k^} is an orthonormal basis for H(E) and 
x = 0Jk for some k E Z. Hence, for such an x, we hx the corresponding 
sequence {<u n } and use the notation x = 0J k . 

Our first step is to rewrite the variational formulation of the Bergman 
metric in the following way: 


Pi(u k ) = 


2 _ sup{|h/(w fc )| 2 : \\h\\ H{E) = 1, h(u) k ) = 0} 


1 


n<P(uk) 

1 

TT(t)'{iOk) 


sup 


sup 


ir 2 K(uJk,Uk) 

ti(uj k ) 2 


E(u k ) 


'■ I HI if (E) — 1' h(u; k ) — 0 


|)'k> 


: IHIh(.e) — 1) h{uik) — 0 


(9) 


Let h E H(E) satisfy IHI-HLE) = 1 and h(ui k ) = 0. Since the functions 

\E(x) \ sin ((j)(x) -(j>(u} n )) 
y/7T(^(u} n )(x - u n ) 

form an orthonormal basis for H(E ), there exists a sequence {c n } E £ 2 (Z) 
such that 

h(x) r k Un (x) r sin (0(x) - <t>{u n )) 

E ( x ) ~t k n E ( x ) ~t k " - w n )e-^( x ) 

= JL = y Cn ( ie ^ 1) AALJ = - L. 

' y/<P(u n ){x - Ul n ) 


As we may assume that the sequence c n is finite, we can differentiate term- 
by-term to get 



(x) 


_j_ ~ ( _ ${%) _A _ f ku n (x) \ 

~^k H \ ~ w n) / n^k n \ E (u k )(x -Un)J ’ 


where c n = c n (e ia (—l) n ). Since k Un (Lo k ) = 0 for every n / k, it holds that 

iV w4e- ^ ^ 




Therefore, the supremum in (J9]) can be rewritten as 

dn. 


Pl(Uk) = 


2 _ (j>'(u k ) 


7G 


sup 


\/${U n )(u k ~ U n 


■■ d k = 0, ^2 \ d n\ 2 = 1 


nEZ 


(10) 
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As this is the dual-formulation of an / 2 -norm. it follows immediately that 


Pi(u k ) 2 = 




7T 




1 


n^k 


4>'(u n )(uj k - u ri 


By Lemma 2.1 we get Theorem 1.5 


Remark 2.2. We observe that, with the notation above, the relation be¬ 
tween the Bergman and the sine distances (see j5j Proposition 8]) gives that 
for any interval / CK with (//-measure equal to 1, i.e., fj (f'(t)dt = 1, then 


J pi(x)dx > 1. 


Therefore, the expected number of real zeros in an interval of bounded <ff- 
measure is bounded below. But these two integrals are in general not equiv¬ 
alent, in the sense that there exist I of (//-measure 1 with fj pi(x)dx = +oo. 
Indeed, this happens with a half-line of (//-measure 1 in the Bessel and Airy 
case (see Remark 3.5). However, we have not been able to find an example 
of a bounded interval where this happens. 


3. Examples 


3.1. The Paley-Wiener space. As we mentioned in the introduction, the 
main example of a de Branges space is the Paley-Wiener space, which cor¬ 
responds to the choice E(z) = e~ mz . We recall that the Paley-Wiener GAF 
is given by 


£ 


sin 7 t[x — n 


7 t(x — n) 


where a n are real i.i.d. standard normal random variables. The phase func¬ 
tion is 4>(x) = 7 tx and trivially 


pi(x) a 4>'(x). 


( 11 ) 


3.2. de Branges spaces with doubling phase function. Here, we con¬ 
sider de Branges spaces with phase function (f>(x) such that the measure 
<f>'(x)dx is locally doubling on M. 

We recall that a measure p is doubling on M if there exists a constant 
C > 0 such that, for all intervals I C M, the inequality p(21) < Cp(I) holds. 
The measure is called locally doubling if it is doubling for all intervals I CK 
satisfying p(I) < 1. 

Locally doubling de Branges spaces are in some sense close to Paley- 
Wiener spaces and are amenable to study. E.g., it is known that these 
spaces correspond (cf. Remark |1.1| ) to so called meromorphic one-component 
model subspaces of the Hardy space, and, in [19j, density-type results for 
sampling and interpolating sequences, analogue to those known to hold for 
Paley-Wiener spaces, were obtained. 

We need the following result, which is Lemma 2.3 in [ T9] . 

Proposition 3.1. Let H(E ) be a de Branges space with phase function </>. 
If the measure <p'(x)dx is locally doubling on M, then there exist constants 
such that 


|(/>(s) - (f(t )| < 1 


4>'(s) ~ 4>'(t). 


( 12 ) 
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As a generalisation of we obtain the following. 

Proposition 3.2. Let H(E ) be a de Branges space with phase function (j). 
If the measure f'(x)dx is locally doubling on M ; then 

Pi(x) ^ 

Before we give the proof of the proposition, we note that by the definition 
of the first intensity function, we immediately get the following corollary. 

Corollary 3.3. Let H(E) be a de Branges space with phase function <f>. If 
the measure <f>\x)dx is locally doubling on M, then, for I cR, the de Branges 
GAF 

F (x ) = ^2a n k U}n (x) 


E[#(Z R (F) n /)] ~ j 4 /(x)dx. 


satisfies 


Proof of proposition \ 3.S\ Define the function 

= / lo sl s (~)l if Im ~ ^ °> 

^ \ log|£*(z)| if ImzCO. 

Let / € H(E ) be such that f(x) = 0 and ||/||_h-(e) = 1- The Bernstein 
inequality (see, e.g., Lemma 2.10 in |19|). says that 

2 ! 

< cfi(x ) 3 / \f(z)\ 2 e~ 2 ^ z) dm(z), 


\f(x)\ 2 

K(x,x) 4>'(x) 


11 (I) 


t D^(x, 1 ) 


where D^(x, 1) is the disc in C, having both diameter and intersection with 
M equal to {t : \(j){x) — cj)(t)\ < 1}. It is not difficult to show that the measure 
(//(Re z )xd 4 ,(x,i) is a Carleson measure for the Hardy spaces H 2 ( C+) and 
H 2 { C_), see [TSjJ Remark 2.12], Therefore 

[ \f{z)\ 2 e~ 2 ^ z) dm(z) < \\f\\ 2 H{E) = 1. 

Jd^x, i) 

Combining the previous two computations, we obtain 

< * w *. 

K(x, x) ~ v K ’ 

Taking into consideration the Bergman metric expression for the first inten¬ 
sity function, it follows immediately that p\(x) < (j)' (x). 

For the opposite inequality, let x, y 6 M be such that (f{x) — 4>(y) = i r, 
and set f(z) = k y (z). Clearly f(x) = 0 and ||/||//(E) = 1- Since 

2 <m 2 \E( X ) \ 2 

1 y{ )l *<ny)\ x - y r 

we deduce, again from the Bergman metric expression for the first intensity, 
that 

\K( x )\ 2 </>'( x ) 


PiW 2 >U 


K{x,x) n(j)'(y)\x - y\ 2 


~ 4>'(xy 


In the last step, we use (12) and the relation <f'(c)(x — y) = 4>(x) — <f(y) given 
by the mean value theorem. ■ 
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3.3. The Airy kernel. The Airy function Ai(z) is the solution of the dif¬ 
ferential equation 


u 

u = zu 


that, for (El, is expressed by the oscillatory integral 


A ’“> = l 


f+oo 3 

e itx+i 1rdx. 


For information about the Airy function one can consult the monograph 
The function E(z) = Ai^z) — iAi(z) =: A(z) — i B(z) is in the Hermite- 
Biehler class and has no real zeros. It therefore yields a de Branges space 
H(E) with reproducing kernel given, for real x y, by 

B(x)A(y) - B(y)A(x ) 


I<(x,y) = 


n(x-y) 

Ai(x)Ai'(y) — Ai'(x)Ai(y) 


(13) 


vr(x - y ) 

The functions Ai(z),Ai / (z) have zeros only on the negative real axis. We 
denote the zeroes of Ai, as is customary, by 

0 > ai > a 2 > .... 

In the following proposition, we summarize some facts about the resulting 
Airy de Branges space. 

Proposition 3.4. Let H(E) be the de Branges space given by the function 

E = Ai'(z) — iAi(z). 

Then, the associated de Branges GAF can be expressed as 

(—l) n+1 Ai(x) 


F{x) = Y J Xn 


n=l 


0r(x - a n ) 


(14) 


where a n are real i.i.d. standard normal variables, and has first intensity 


1 


Pi{x) = -o o 


Ai(x)Ai'(x) 


Ai 4 (x) 


x 


7 r 2 \3xAi(x) 2 — Ai'(x) 2 4(xAi(x) 2 — Ai'(x) 2 ) 2 3 

Moreover, if we denote by cf the phase function of E, the following holds: 

(i) The measure 4>'(x)dx is not locally doubling on M. 

(ii) The function <//(x) satisfies the estimates 


(15) 


fi'(x) 


+ 0(x 5 / 2 ) 


2a; 3 / 2 

i(sinR||xr-|) + 


as x 


+oo, 


cos 2 (||a:| 3 / 2 — 
1*1 


-1 


(1SI-I —oo. 


(iii) The first intensity function satisfies the estimates 

. . f ik + 0(a; _3/2 ) as x -A +oo, 

- i /w ^ 

- \ V as x -A —oo. 

V 7T V 3 

One can compare this estimate with the corresponding one for the the 
determinantal case ( (281 (1.17)]) where the first intensity behaves as 7r~ 2 ^/jx| 
when x —> —oo and decreases faster than exponentially when x —> +oo. 
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Figure 1. First corre¬ 
lation function, pi(x). 



Figure 2. Derivative 
of the phase function, 
(j)'(x). 


Remark 3.5. If we compare cp 1 with p \, we observe the following. For x < 0, 
the first intensity function behaves, in average, like the oscillating function 
<j>'{x). Indeed, for I C (—oo,0), it holds that 


J (p'(x)dx = 1 =i> J pi(x)dx — 1 

To see this, simply use the well-known estimate (see 

2/3 

a„ = - ( I ( 1 + o 


formula (2.52)]), 


(^) M 1 *)' 


to see that 

r a k r(k+ l ) 2/3 

/ Pi(x)dx ~ / y/xdx ~ 1. 

Jak+i Jk 2 / 3 

For x > 0 the situation is quite different. Indeed, the analogy that can 
be seen when taking averages breaks down for x > 0, since, unlike in the 
determinantal case, the expected number of zeros on (0, oo) is unbounded 
while it has finite <//- measure . Indeed, the expected number of zeroes on the 


interval (0, IV) grows as log N (see Remark 2.2). 


Proof. Since Ai and Ai r do not have any zeroes in common, it follows that 
E is real exactly on the points, a n , and so, these form a sequence of points 
which are exactly at (//-distance 7r from each other. Since \Ai(x)/E(x)\ 2 is 
not integrable (that it is not integrable on (0, oo) frollows from [UJ items 
10.4.59 and 10.4.61, p.448] or [32] formulas (2.44) and (2.45)]), by Lemma 


1.2 the normalized reproducing kernels {k(x, On)}?ieN yield an orthonormal 


basis of H(E). 

To obtain an explicit formula for the Airy de Branges GAF, we note that 
from the relation Ai"(x) = xAi(x), we get 

K(x,x) A(x)B'{x) — A'(x)B(x) 


<j)'{x ) = 7T 


\E(x)\ 2 A 2 (x) + B 2 (x) 

Ai'(.T) 2 — xAi(x) 2 


(16) 


Ai'(x) 2 + Ai(x) 2 

Taking into account that Ai(a n ) = 0 and Ai(x) is real-valued on M, it is now 
straight-forward to see that in the Airy case, formula ([!]) yields (14). 
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To compute p\ for the Airy de Branges GAF, it is practical to use the 
following representation for the reproducing kernel 


I\(x,y) = - 
7 r 


r »+00 


Ai(x + t)Ai(y + t)dt, 


that follows from |32l formula (3.52)] and the fact that both Ai(x) and Ai^x) 
decay as x —> +oo. Combined with the Edelman-Kostlan formula Q and 
formulas for primitives of products involving the Airy function, [32] Chapter 


3], it yields (15). We remark that one can get a similar formula for the first 
intensity function for any de Branges space considering the related Sturm- 
Liouville (or Schrodinger) equation (see [23.1T8]). 

To see that the Airy space is not locally doubling, we verify that ( |12| ) does 
not hold. To this end, denote the (negative) zeros of Ai 7 (x) by b k . Then (see 
[ 32] p. 15]), it suffices to observe that 

\</>(bk) - 4>{a k )| = and (j)'(a k ) = &(b k ) = — ~ k 2/3 . 

2 7T 7T 

More generally, one may translate asymptotic estimates for the Airy func¬ 
tion and its derivative in combination with (16) and (15) to check that (ii) 
and (iii) holds, respectively. More specifically, in both cases, one uses the 
estimates [TJ 10.4.59 and 10.4.61] to obtain the asymptotics asx-> +oo and 
the estimates [1, 10.4.60 and 10.4.62] to obtain the asymptotics as x —> — oo. 


3.4. The Bessel kernel. For v > —1/2, the Bessel function of order v, 
denoted by J u , is the solution to the differential equation 

x 2 y" + xy' + {x 2 - v 2 )y = 0, (17) 


that can be expressed, using Taylor series, as 

00 l i\k(z\2k+v 

j ( z ) = v ^ [<2) _ 

^ k\T(k + v + 1)' 

k=0 

One can define a de Branges space using the Hermite-Biehler function (see 

m) 


E u (z) = A u (z ) - i B u (z) 

Observe that with the term z~ v ^ 2 in the definition, the function E u (z) is 
entire. The reproducing kernel for the de Branges space H(E U ) is given by 


K(x, y) 


{xy) 


„/ 2 Vy J 'v(Vy) J v(Vx) - VxJ' v {Vx)j u {^y) 

tt(x - y) 


x,i/eR. 


(18) 

In the determinantal case, it is more common to define the kernel without 
the extra term ( xy )~ u / 2 , but this version appears also in the literature. See, 


e.g., |H]. 

The functions J u and J' u only have zeroes on the positive half of the real 
line, and we denote the zeroes of J u by 


0 < ju, l < ju, 2 < • ■ • • 
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We summarize some facts about the Bessel de Branges GAF in the following 
proposition. 

Proposition 3.6. For v > —1/2, let H(E U ) be the de Branges space given 
by the function E v defined above. Then the associated de Branges GAF can 
be expressed by 

F(x) = Jl X -» Y. “» { ~ i)kj ‘" kJ fi VS) , (19) 

V vr x — 7 , 

k =l 

where a n are real i.i.d. standard normal variables 

Moreover, if we denote by the phase function of E u , the following holds: 

(i) The measure f>'(x)dx is not locally doubling on M. 

(ii) The function 4>'{x) satisfies the estimates 

+ 0 (x -1 / 2 ) as x +oo, 


as x 


—oo, 


it(\ I 2 (y/x sin t+C 2 cos f) 2 + cos 2 t—2c± sin 2 t 

h_v ! + 0(i - 2) 

where t = \Jx — ^ /i = An 2 , c\ = (A ' am ^ c 2 = 

(iii) The first intensity function satisfies the estimates 

( i 


Ma H 3T+o(i*iw » 


as x +oo, 
x —> —oo. 


Observe that the formula for (fi above says that asr-> +oo, we have that 
2 cj/(x) oscillates between l/(x — 2 ci) and the constant 1/(1 + c|). 

As in the Airy case, we can interpret the first intensity function as giving 
an “average” of the (//-function, since, for k > 1 , we have 


rE,k +1 

•A? 2 ! 


pi(x)dx 




4 2 

Jv.k 


dp (x)d'x = n. 


And, as in the Airy case, the analogy breaks for on the other half line: 

ril, i fil, i 

n = (p'(x)dx < / p\(x)dx = oo, 

J— OO J— OO 

see Remark 12.21 

Also, it is not difficult to estimate the first intensity for the the determi- 
nantal process given by the kernel K(x, y)(xy) u ^ 2 


K(x,x)x u ck ^ 1/2 ^1 + Q(x 1/2 )^ 


for t = \fx — yy — f, and we get a bigger amount of points (in expectation) 
than in the GAF case. 

Proof. Since E v is real exactly on the points, j v ,k-, these are necessarily at 
(//-distance 7 r from each other. As in the Airy case, by using the well-known 
asymptotic estimates for the Bessel function and its derivative m items 
9.2.5 and 9.2.11, p. 364]), one can see that \B u (x)/E u (x)\ 2 is not integrable 
and therefore, by Lemma 


1.2 


m 


that {K(x,j 2 k )}k >i is an orthogonal basis i 

H\E). 

A straight forward calculation yields that the Bessel de Branges GAF may 
be given by the formula (19). 
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Figure 3. First corre¬ 
lation function p\ (x ), 
for v = 


Figure 4. Derivative 
of the phase function 
for v = -. 


By using Bessel’s equation ( fl7| ), one can show that 

B' v {x)A v {x) - B v (x)A' v {x) 




1 


A 2 {x)+B 2 (x) 

JUV ^) 2 + (i - #) MVzf 


( 20 ) 


2 xJ' v (y/x) 2 + J„(y/x) 2 

From this, we can observe, as in the Airy case, that the de Branges spaces 
given by the Bessel functions are not locally doubling. Indeed, since the 
zeros of J u and J' interlace (i.e., j v ^ < j' v k+l < j v ,k+ 1, where we let j' vk 
denote the zeroes of J'), one can make the following argument. If H(E U ) 
was locally doubling, then, by property (12), it would follow that 


ML) * ML) 


From the expression of i ft, we get 


<Aih) = 

z Ju,k 


and 


ML) 


1 


V \J V ,k> 


2 Jvtil 




Since J v (x) satishes Bessel’s equation of order u, i.e., equation ©, we get 
a contradiction from the limit 

Wujk) 1 

J 7 ./ x = ^ 2 ~ - 1 fc->oo. 

Wu,k) t,k 

Compared to the Airy case, it is much more difficult to get asymptotic 
estimates for the phase function, its derivative, and the first intensity func¬ 
tion. 

Using formula ( |20[ ) and an accurate estimate for J v (x) and its derivative 
J' v {x) (see, e.g., [U items 9.2.5 and 9.2.11, p. 364]), we get the estimate for 
4>'{x) as x —> +oo. To estimate what happens when x —> — oo, we combine 


(20) with asymptotic estimates for the modified Bessel functions of the first 


kind I v {z) and its derivative (see, e.g., jU items 9.7.1 and 9.7.3, p. 377-378]). 
Here, we also need the relation 


/ \ UL7V / \ 
J v [iz) = e 2 I u (z), 


which is valid for —7r < arg(z) < and so, in particular for z = y/x for 
x > 0. 
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To obtain an explicit expression for pi(x), in principle it is possible to 
proceed in the same way as for the Airy function, using the reproducing 
property [3U P- 295] 

2ttK(u, v) = (uv)~ U//2 f J u (Vut)Jv(Vvt)dt 

Jo 

and the Edelman-Kostlan formula. However, this requires far more work 
than in the Airy case. Instead, it is possible to get the estimate as a; -A +oo 
using direct estimates on ft, while to get an estimate as x —> — oo, one can 
combine direct estimates with Lemma 11.71 ■ 


3.5. Orthogonal Polynomials. The de Branges spaces above all belong to 
a more general class of spaces which can be built in connection with singular 
or regular Sturm-Liouville problems. More precisely, the Airy equation on 
[0, Too) and the Bessel equation in (0,1] are examples of singular Sturm- 
Liouville problem in the limit-point and in the circle-point case. For more 
details see [231 EM- A discrete version of these problem is to consider spaces 
of orthogonal polynomials. 

Let /i be a finite positive Borel measure p on M whose support is not a 
finite number of points and which has finite moments of all orders, i.e., 

/ \x\ n dp < oo, Vn G No- 

Jr 


Denote by p n (x) = r y n x n T • • •, 7 n > 0, a set of polynomials orthonormal 
with respect to the measure p, i.e., satisfying 


r+OO 


pn{x)p m {x)dp(x) = 5 n 


One can define the Hermite-Biehler function E n (z ) = y7r 7 " 1 (p n _\(z) — 

ipn(z)), and there is a unitary equivalence between the de Branges space 
H(E n ) and the space of orthonormal polynomials with respect to p. The 
associated de Branges GAF has the same distribution as the random poly¬ 
nomials 

71—1 

Pn(x) = ^ a kPk(x) (21) 

k =0 

where a*, are real i.i.d. real standard normal random variables. This is clear 
because the covariance kernel of this gaussian proces is given by K n (x,y ) = 
TJk=oPk{x)Pk{y)- 

The problem of estimating asymptotically the expected number of real 
zeros has been studied in some cases, see [33] for the Legendre case. 

In principle, by identifying the phase, it would be possible to compute the 
expected number of real zeros, but this seems far from being trivial. 

One case which can be studied explicitly is the orthonormal polynomials 
with respect to the measure d^ = (■x 2 T a 2 )~ n dx, where a > 0 and n > 2. 
In this case, one can take E{z) = (z T ia) n . If 4>(x) is a primitive of the 
Poisson kernel, i.e. 4>'(x) = j3 ° a2 , a straight-forward computation yields 
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50 = — 2(0') 2 . From this it follows that a phase function for H(E) is n0(x) 
and 


Pi{x) = ~ 

7r 


1 n 2 - 1 


4>'{x). 


Consequently, the expected number of real zeros of (21) is, in this case, 


E[#(z K (P n )ni)] = / Pl {x) = 


n 2 — 1 


4. Proofs of Theorems PI AND P 

We begin with a lemma on the Schwarzian derivative. It is one of four 
lemmas from which Theorem 11.91 follows. 


Lemma 4.1. Suppose that, for all i£l, 

S(j) i(x) + 20 / 1 (x) 2 = 50 2 (x) + 20' 2 (x) 2 . 


( 22 ) 


50i(0i 1 (t)) = - — f ' 


If we put if = 02 o (f > 1 1 , then it holds that 

1 = ijj'{t) 2 + ^Sip(t). 

Proof. The Chain Rule for the Schwarzian derivative is 

S(f°g) = (Sf o g){g') 2 + Sg, 

whence 

s(^m 

Combined with the usual chain rule, this yields 

2 ( f , 1 Uf\t)) 2 + 50i(00 1 (t)) =-^- SWiW) 

Yl ^i w/ (cf^nty (^nt) 2 

Applying the hypothesis to the left-hand side above, with x = 0(j” 1 (t), and 
rearranging the terms, this becomes 

2 = 2</>’ 2 (<l>i 1 (t)) 2 (<f>i 1 )'(t) 2 + 50 2 (0r 1 (*))(0 f 1 )'^) 2 + 5(00 1 )(t). 

Now, using both chain rules on the right-hand side, we get 
2 = 2(^2 o 00 1 )'(f) 2 + 5(02 o 0r 1 )(O, 

whence the result follows. ■ 


Introducing the change of variables x = log if 1 and y 


(22) into a second order differential equation. Indeed, 


f x' =y, 

\ y' = 2 — 2e 2x + \y 2 . 


x', we can turn 


(23) 


To obtain a linear differential equation, in terms of differentiation with 
respect to x, we observe that, by the chain rule, 


dy 

dx 


dy dt 
dt dx 


2 - 2e 2x + \y 2 
V 

< 2 - 2^1 + 0 ,. 
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Figure 5. Phase portrait of (23). 


It now follows, by the further change of variables z = y 2 , that 


dz 

dx 


= z + 2(2 — 2e : 


2x\ 


This linear differential equation has the general solution 

z(x) = Ce x — 4(1 + e 2x ), 

whence, by substituting back, it follows that if (x(t),y(t)) is a solution of 
(23), then 

y(t) 2 = Ce x ®-4(1 +e 2 *®). 

Therefore, the system has periodic orbits, and, by taking different values of 
C > 8, one gets the phase portrait in Figure [5] 

The next result shows that the function (t) is 7r-periodic. 


Lemma 4.2. The solutions to the system (23) have period n. 


Proof. Let (x, y) be a solution to (23). As x' = y, it is clear that it is enough 
to show that x(t) has period n. For C > 8, let x± = x±(C) be the two 
numbers such that Ce x — 4(1 + e 2x ) = 0, i.e., the points where y = 0 in the 
phase portrait. 

Let T be the period of the solution (x,y). If we let x(0) = x_, then, by 
symmetry, from time 0 to time T/2 the function x(t) goes from x_ to x+. 
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Therefore, by the inverse change of variable x(t) = s, we get 

2 Jo J x _ ^Ce° - 4(1 + e 2s ) 


rr+ 


dr 


n/r+ TyJCr - 4(1 + r 2 ) 


arcsm 


Cr — 8 


1 r+ 


7r 


VC^8VJ 1/r+ 2- 

Here, we used, in the first line, that 

x'(t) = y(t) = \J Ce x W — 4(1 + e 2x JI) = \/(7e s — 4(1 + e 2s ), 

and, in the second line, the change of variables e s = r. Since x± corresponds 
to y = 0, it follows that 

1 (C VC 2 -8 2 \ 

r± = 2(T*—i—]■ 

and x_|_ + X- = logr+ + logr_ = log r+r_ = 0. The last step of the 
computation now follows. ■ 


Lemma 4.3. Let (x,y) be a solution of the system (23) such that y(0) = 0. 

If 


U(s) = [*e x ®dt, 

Jo 


then, for all k G Z, we have 


Proof. By essentially the same computations as in the previous proof, we get 


U[k~) = k—. 


U{r r/2) = 


rvr/2 


^dt = 


r r + 


dr 


<l/r+ \/Cr — 4(1 + r 2 ) 
f r+ du 7r 

Jl/r + uy/Cu — 4(1 + u 2 ) 2 

Here, the only difference is the change of variables u = l/r in the second 
line. 

With this, it is enough to show that, for 0 < t < n, 


Indeed, defining 


x(t) = x(n — t ). 


(x(t),y(t)) = (x(?r - t), -x'(tt - t)), 


we see that (x,y) is a solution of (23) with initial condition (x(0),y(0)) = 
(xo,0). Therefore, it follows by uniqueness that 

(x(t),y(t)) = (x(tt - t ), -x'(tt - t)), 


whence the desired conclusion follows. 
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The above result shows that if (02 o (f 1 1 ) /, (0) = 0, which is tantamount to 
y( 0) = 0, then the function 


/ (</>2 ° 4 >i 
Jo 


is 7r-periodic. 

Let a € M be such that (02 0 0= 0. Indeed, such an a exists 
because (log 00 0 /; /0 7 ) for 0 = 02O00 1 is a solution of the system (23). Now, 
we observe that the solution (x,y) to (23), where x0) = log(02o00 Yi^S-a), 
satisfies x'(0) = 0. Hence, by Lemma 4.3 above, 

(02 ° 0c 1 ) (y + ~ 0 0r 1 )(«) = -0“- 

So, for Vk = 05~ 1 ( z y + a), we get 

02 04c) - 02 04)) = 01040 - 0104))- 
Therefore, we have that 

02040 = 01040 + 0, 

where {u/0 are the points such that 0i(u/0 = a (mod Ip, and /3 £ 1 is a 
constant. 


The proof of Theorem 1.9 is now a consequence of the following modi¬ 
fication of a result of de Branges [9j Theorem 24], which is well-known to 
specialists. 


Lemma 4.4. Let E\ and E 2 be two Hermite-Biehler functions without real 
zeroes, with phase functions 0i and 02, respectively. Suppose that there exist 
a, /3 G M so that 

02 0 ) = 010 ) + 0 

for all t E R. such that 0i0) = a(rnod^). Then, there exists a non van¬ 
ishing real entire function S such that F{z) 1 —> S(z)F(z) is an isometric 
isomorphism from Fd{Ei) onto H{E 2 ). 

Proof Define the Hermite Biehler functions E\ and E 2 by 

E^z) = e-^-^E^z) , and E 2 (z) = e ia E 2 (z). 

If 0i, and 02 denote their corresponding phase functions, then 

010 ) = 02 0 ), 

whenever t is such that 02 0) = 0 (mod |). Therefore, by Theorem 24 in (9], 
there exists a real entire function S(z) such that 

F{z) ha S(z)F(z) 

is an isometric map between the spaces H(E\) and L0.&0. Since H{E\) = 
H(Ei) and H{E 2 ) = H{E< 2 ), in the sense of Hilbert spaces, it follows that S 
induces an isometry between the original spaces H(E\) and H{E 2 ). A priori, 
this function may have real zeroes. However, since neither E\ nor E 2 have 
real zeros, the function S never vanishes. Indeed, the function S satisfies the 
identity 


K 2 (z,w) = S(z)Ki(z,w)S(w) 
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and Kj(x, x) = ^'•(x)|£' j (V)| 2 for j = 
that, for every iGl, 


\S(x)\ 


K 2 (x,x) 

K\(x,x) 


1, 2 and i£l. From this it follows 
f 2 {x)\E 2 (x)\ 2 


From this we conclude that S is an isometric isomorphism between the spaces 
H{E 1 ) and H(E 2 ). ■ 


Now we can easily deduce Theorem 1.8 from Theorem |1 .9 


Proof. (Theorem 1.8) Let F{x ) be a de Branges GAF defined by the repro¬ 
ducing kernel I\i(z,w) of the space H(E\) and let G{x) be the de Branges 
GAF defined by the reproducing kernel K 2 (z,w) of the space Ff{E 2 ). It fol¬ 
lows from Theorem |1.9| that S(x)F(x) and G(x) have the same covariance 
kernel and therefore, as they are gaussian processes, S(x)F(x) and G(x) 
have the same distribution, but S(x) does not vanish, so we get the result. 


References 

1. Milton Abramowitz and Irene A. Stegun, Handbook of mathematical functions with 
formulas, graphs, and mathematical tables, National Bureau of Standards Applied 
Mathematics Series, vol. 55, For sale by the Superintendent of Documents, U.S. Gov¬ 
ernment Printing Office, Washington, D.C., 1964. 

2. Robert J. Adler and Jonathan E. Taylor, Random fields and geometry, Springer Mono¬ 
graphs in Mathematics, Springer, New York, 2007. 

3. Lars V. Alilfors, Cross-ratios and Schwarzian derivatives in R’ 1 , Complex analysis, 
Birkhauser, Basel, 1988, pp. 1-15. 

4. Jorge Antezana, Jeremiah Buckley, Jordi Marzo, and Jan-Fredrik Olsen, Gap proba¬ 
bilities for the cardinal sine, J. Math. Anal. Appl. 396 (2012), no. 2, 466-472. 

5. N. Arcozzi, R. Rochberg, E. Sawyer, and B. D. Wick, Distance functions for repro¬ 
ducing kernel Hilbert spaces, Function spaces in modern analysis, Contemp. Math., 
vol. 547, Amer. Math. Soc., Providence, RI, 2011, pp. 25-53. 

6. Stefan Bergman, The kernel function and conformal mapping, revised ed., American 
Mathematical Society, Providence, R.I., 1970, Mathematical Surveys, No. V. 

7. Gordon Blower, Random matrices: high dimensional phenomena, London Mathemat¬ 
ical Society Lecture Note Series, vol. 367, Cambridge University Press, Cambridge, 
2009. 

8. Eugenio Calabi, Isometric imbedding of complex manifolds, Ann. of Math. (2) 58 
(1953), 1-23. 

9. Louis de Branges, Hilbert spaces of entire functions, Prentice-Hall, Inc., Englewood 
Cliffs, N.J., 1968. 

10. Alan Edelman and Eric Kostlan, How many zeros of a random polynomial are real?, 
Bull. Amer. Math. Soc. (N.S.) 32 (1995), no. 1, 1-37. 

11. Naomi D. Feldheim, Zeroes of Gaussian analytic functions with translation-invariant 
distribution, Israel J. Math. 195 (2013), no. 1, 317-345. 

12. Victor Havin and Javad Mashreghi, Admissible majorants for model subspaces of H 2 . 
I. Slow winding of the generating inner function, Canad. J. Math. 55 (2003), no. 6, 
1231 -1263. 

13. J. Ben Hough, Manjunath Krishnapur, Yuval Peres, and Balint Virag, Zeros of Gauss¬ 
ian analytic functions and determinantal point processes, University Lecture Series, 
vol. 51, American Mathematical Society, Providence, RI, 2009. 

14. Marek Jarnicki and Peter Pflug, Invariant distances and metrics in complex analysis, 
extended ed., de Gruyter Expositions in Mathematics, vol. 9, Walter de Gruyter 
GmbH & Co. KG, Berlin, 2013. 






24 


JORGE ANTEZANA, JORDI MARZO, AND JAN-FREDRIK OLSEN 


15. M. Kac, On the average number of real roots of a random algebraic equation, Bull. 
Amer. Math. Soc. 49 (1943), 314-320. 

16. H. Blaine Lawson, Jr., The Riemannian geometry of holomorphic curves, Bol. Soc. 
Brasil. Mat. 2 (1971), no. 1, 45-62. 

17. D. S. Lubinsky, Universality limits at the hard edge of the spectrum for measures with 
compact support, Int. Math. Res. Not. IMRN (2008), Art. ID rnn 099, 39. 

18. N. Makarov and A. Poltoratski, Meromorphic inner functions, Toeplitz kernels and the 
uncertainty principle, Perspectives in analysis, Math. Phys. Stud., vol. 27, Springer, 
Berlin, 2005, pp. 185-252. 

19. Jordi Marzo, Shahaf Nitzan, and Jan-Fredrik Olsen, Sampling and interpolation in de 
Branges spaces with doubling phase, J. Anal. Math. 117 (2012), 365-395. 

20. Madan Lai Mehta, Random matrices, second ed., Academic Press, Inc., Boston, MA, 
1991. 

21. Nikolai Nikolski, Distance formulae and invariant subspaces, with an application to 
localization of zeros of the Riemann (-function, Ann. Inst. Fourier (Grenoble) 45 
(1995), no. 1, 143-159. 

22. Brad Osgood, Old and new on the Schwarzian derivative, Quasiconformal mappings 
and analysis (Ann Arbor, MI, 1995), Springer, New York, 1998, pp. 275-308. 

23. Christian Rending, Schrodinger operators and de Branges spaces, J. Funct. Anal. 196 
(2002), no. 2, 323-394. 

24. S. O. Rice, The Distribution of the Maxima of a Random Curve, Amer. J. Math. 61 
(1939), no. 2, 409-416. 

25. _, Mathematical analysis of random noise, Bell System Tech. J. 24 (1945), 46- 

156. 

26. M. Sodin, Zeros of Gaussian analytic functions, Math. Res. Lett. 7 (2000), no. 4, 
371-381. 

27. A. Soshnikov, Determinantal random point fields, Uspekhi Mat. Nauk 55 (2000), 
no. 5(335), 107-160. 

28. Alexander B. Soshnikov, Gaussian fluctuation for the number of particles in Airy, 
Bessel, sine, and other determinantal random point fields, J. Statist. Phys. 100 (2000), 
no. 3-4, 491-522. 

29. Terence Tao and Van Vu, Local universality of zeroes of random polynomials, 
arXiv: 1307.4357 [math.PR], 2013. 

30. Craig A. Tracy and Harold Widom, Level-spacing distributions and the Airy kernel, 
Comm. Math. Phys. 159 (1994), no. 1, 151 174. 

31. _, Level spacing distributions and the Bessel kernel, Comm. Math. Phys. 161 

(1994), no. 2, 289-309. 

32. Olivier Vallee and Manuel Soares, Airy functions and applications to physics, second 
ed., Imperial College Press, London, 2010. 

33. J. Ernest Wilkins, Jr., The expected value of the number of real zeros of a random sum 
of Legendre polynomials, Proc. Amer. Math. Soc. 125 (1997), no. 5, 1531-1536. 

Universidad Nacional de La Plata, Departamento de Matematica, Esq. 50 

y 115 s/n, Facultad de Ciencias Exactas de La Plata (1900), Buenos Aires, 

Argentina 

E-mail address : antezana@mate.unlp.edu.ar 

Departament de matematica aplicada i analisi, Universitat de Barcelona, 

Gran Via 585, 08007, Barcelona, Spain 

E-mail address : jmarzo@ub.edu 

Centre for Mathematical Sciences, Lund University, P.O. Box 118, SE-221 

00 Lund, Sweden 

E-mail address', janfreol@maths.lth.se 





